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Abstract 

> 

We study the measurement of the quark fragmentation functions for spin-1/2 

• baryon production (A and A, in particular) in unpolarized e + e~ annihilation. 
The spin-dependent fragmentation functions g\{z) and h\{z) can be probed 
in the process as a result of quark-antiquark spin correlation and the weak 

<3\ ' decay of the baryons. The relevant cross section is expressed as a product 

• of the two-jet cross-section, the fragmentation functions, and the differential 
width of the hyperon decay. 
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I. INTRODUTION 



Due to confinement, the basic building blocks of QCD, quarks and gluons, cannot emerge 
as asymptotic states of the theory. Rather the high energy quarks and gluons in a hard 
scattering process show up in a detector as jets of hadrons, and their axes register the 
momentum directions of the initial partons. The process of converting a colored parton 
into a shower of hadrons is called parton fragmentation, in which it has long been known 
that the soft QCD physics dominates. As a result, analyzing fragmentation directly from 
first principles is very difficult. Up to now the fragmentation process is mainly modelled 
with a set of empirical rules and a large Monte Carlo code [|TJ. Given the paucity of our 
knowledge, experimental information on various aspects of fragmentation is quite valuable 
for understanding the essential physics in the underlying process. 

Semi-inclusive information about fragmentation is contained in fragmentation functions, 
which measure the probability for a quark to fragment to a specific hadron with a fixed 
longitudinal momentum. All the fragmentation functions experimentally measured so far 
are spin-independent. In Refs. [@-@], the novel spin-dependent fragmentation functions gi(z) 
and hi(z) are introduced for spin- 1/2 baryon production. Both g~i(z) and hi(z) are of leading 
twist and can produce significant effects in high-energy processes. The physics of these 
fragmentation functions is quite simple. The g\{z) fragmentation function represents the 
probability density to find a spin-1/2 baryon with its polarization in the direction of the 
longitudinal polarization of the original quark. Likewise, hi(z) represents the probability 
density to find a spin-1/2 baryon with its polarization in the direction of the transverse 
polarization of the original quark. Here, longitudinal or transverse is defined relative to the 
momentum of the quark. Undoubtedly, experimental data on these fragmentation functions 
contain important information about spin-transfer during quark fragmentation. 

To measure the spin-dependent fragmentation functions, one normally needs to produce a 
polarized quark jet. This is possible with polarized e + e~ collisions or polarized deep-inelastic 
scattering. On the other hand, in an unpolarized e + e~ annihilation, the polarizations of 
the created quark and antiquark pair are correlated due to the chiral invariance of the 
electroweak interactions. Therefore if one can find a way to measure the spin correlation 
of the produced hadrons, one can obtain information on the spin-dependent fragmentation 
functions. This observation forms the foundation for our work. 

In this paper, we study the measurement of the spin-dependent fragmentation functions 
in unpolarized e + e~ annihilation. For definiteness, we consider correlated production of 
A and A in quark and antiquark jets respectively. The A and A are observed by their 
weak decay to nucleon-pion and anti-nucleon-pion respectively. The weak decay allows one 
to measure the complete spin-density matrix of the parent baryon. The spin- dependent 
information is extracted from the angular correlation of the decay products of A and A. 
Our result, of course, is valid for production of any spin-1/2 baryon which decays weakly to 
another spin-1/2 baryon and a pion. 

This paper is organized as follows. In §11, we write down the differential cross section as 
a product of the density matrices for the jet production, the quark fragmentation, and the 
baryon decay. We then present the helicity formalism on which the subsequent calculation 
is based. In §111, we evaluate the jet production density matrix. In §IV, the quark frag- 
mentation density matrix is defined and expressed in terms of the fragmentation functions 
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fi(z), gi(z), and h\(z). In §V, A and A decay density matrices are constructed. And finally 
in §VI, we put everything together to form the final cross section, and discuss its physical 
significance. We summarize our results in §VII. The reader who is interested primarily in 
our results and their physical interpretation might omit §111, §IV, and §V and skip directly 
to the formulas and discussions in 5 VI and 8 VII. 



II. SPIN DENSITY-MATRIX FORMALISM IN HELICITY BASIS 

The calculation of the cross section for e~e + — > qq — > AXAX — > p-K~ Xp-K + X can be 
carried out by following the standard procedure of the operator production expansion, or 
equivalently, by using the collinear expansion and "cut diagram" technique However, 
since we are interested in a result at the leading twist, there is a more physical approach 
which employs the language of the parton model and follows the process step by step. First, 
we calculate the quark-antiquark production; then, we deal with quark and antiquark frag- 
mentations; next, we include the A and A decay; and finally, we assemble the qq production, 
fragmentation and decay processes together to obtain the experimentally observable particle 
distributions. The only subtlety in the case of a spin-dependent calculation is that all the 
intermediate quantities must be in a form of spin-density matrices. Then the cross section 
is expressed as a trace of product of spin-density matrices for separate subprocesses. In our 
case we have, 
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where h(h) and H(H) are indices labelling spin states of quark (antiquark) and hyperon 
(antihyperon), respectively. By convention, repeated indices are summed over. The bulk of 
this paper is devoted to defining and calculating these spin-density matrices. 

One can choose any spin-density matrix formalism to perform the calculation. The 
simplest is perhaps the one in which all the spin indices are just the ordinary Dirac indices. 
However, the drawback of this approach is that all the density matrices are 4x4 and they 
do not have a clear physical interpretation. As is often the case, the physics is much clearer 
in the helicity basis 0. Throughout this paper, we shall use this formalism, and thus h(h') 
and H(H') in Eq. (|XJ) are to be interpreted as helicity indices with values ±1/2 (or ± for 
short). 

In the helicity formalism, one first choose the helicity basis, u(h) for fermion and v(h) for 
antifermion. Then a general polarization state for fermion (antifermion) can be expressed 
as, 

U = au(+) + bu(-); V = cv(+) + dv(-) (2) 
or simply in terms of the two-component spinors, 

U=(l); V=( c d ). (3) 
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Any processes and subprocesses with fermions and/or antifermions as external particles can 
be calculated as a density matrix with pairs of helicity indices. Each pair corresponds to 
one particle, with one index representing the spin state of the particle in the amplitude and 
the other in the conjugate amplitude. Obviously the spin-dependent probability or cross 
section can be obtained by contracting these indices with the two-component wave functions 
in Eq. @. 

The helicity formalism treats the polarization of each fermion (antifermion) indepen- 
dently, since the two-component wave function is expanded in basis states which are different 
for each particle. In particular, the spin quantization axis is the direction of the particle's 
momentum. This is rather convenient because one can quickly write down a two-component 
wave function once one knows the relative orientation between the particle's momentum and 
its spin vector. In fact, it will be useful for us to define a coordinate system for each particle 
with its momentum as the z-axis. The directions of the x and y axes are defined as the 
directions of the spin vector S 71 = uj^j^u associated with the following spinors, 

u{x) = -j= («(+) + «(-)) , 

u(y) = -^(u(+) + iu(-)) . (4) 

This way one is assured that the x and y axes are defined in the same way for all particles 
and antiparticles. The choice of helicity basis for each particle fixes the relative orientation 
of the different coordinate systems. For fermions, the two-component spin wave function 
can be chosen universally as, 



a \ ( cos g 



bj V sin §e^ 



(5) 



where 9 and <ft are polar and azimuthal angles of the spin vector relative to the momentum 
and the a;-axis. 

Another important advantage of the helicity formalism is that an antifermion can be 
treated exactly like a fermion. This should be the case because the definition of fermion and 
antifermion is itself arbitrary due to charge conjugation symmetry. However, in ordinary 
calculations, some aspects of the symmetry are not obvious because antifermions are treated 
as holes in the negative energy fermion sea. For instance, the spinor v is associated with 
antifermion creation, whereas u is with fermion annihilation. In other words in a Feynman 
diagram the momentum flow for antifermions is always against the fermion number flow, 
and the Dirac algebra follows the latter flow. In the helicity formalism, if one were to 
follow the fermion number flow in ordering the helicity indices, one would find that all the 
density matrices would be expressed naturally in terms of the transpose of Pauli matrices. 
Furthermore, if one were to use the charge conjugation relation v = Cu T to define the spinor 
for the antifermion, then the two-component wave function would be the complex conjugate 
of that in Eq. (§), 



C\ ( COS ~ 



d J V sin fe-** 



(6) 



Note again that this wave function is associated with the antifermion creation. All these 
observations point to a simple scheme for dealing with antifermion spin calculations: Order 
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the spin indices in terms of the momentum flow rather than the fermion number flow. Then 
the spin-density matrices are all naturally expressed in terms of Pauli matrices and the 
two-component wave functions of antifermions are exactly the same as those of fermions. 
The wave function in Eq. (||) now appears in a transposed form and is naturally associated 
with an antiparticle being created. One final note: the x and y axes for the antiparticle are 
defined by the spin vector S^ L = —vj^j^v associated with the spinors, 

v(x) = -j={v{+)+v{-)) , 

v(y) = ^= («(+) - iv(-)) , (7) 

where the minus sign follows from Eq. (|6|). 

The following sections provide concrete examples for illustrating the helicity formalism 
in detail. 



III. QUARK- ANTIQUARK PRODUCTION DENSITY MATRIX 



Our calculation begins with quark- ant iquark production via photons (7*) and Z°s. Here 
we work in the rest frame of 7* or Z°. Event by event, we choose the z axis in the direction of 
the quark jet, so the anti-quark momentum is in the — z direction. [The identification of the 
quark jet can be made through study of jet ensembles. However, in the present calculation 
we implicitly assume that the A is a fragment of the quark. The possibility of q — > AX can 
be included in the final result in a straightforward way, but the A's resulting from q —>■ AX 
will dominantly occur at low z (the longitudinal momentum fraction of the quark carried 
by A) and usually will not satisfy jet isolation cuts. Henceforth we ignore the possibility of 
q — > AX fragmentation.] By convention, the y axis is defined by y = k x z, where k is the 
direction of the electron beam. In this coordinate system, the 4-momenta of electron and 
positron are 



fc" = (E,-E smQ,0,EcosG) , k' 1 * = (E,EsinQ, 0, —Ecos O), 



(8) 



respectively, where O is the polar angle of the electron. 

In this section, we seek an expression for the quark- ant iquark production density matrix 
(da(e~e + — > qq) / dQ q ) hh , rr, in the helicity basis, where the ordering of the indices follows 
the momentum flow of quark and antiquark, and is shown explicitly in Fig. 1. While the 
diagonal elements of the density matrix are unique, the off-diagonal ones depend on the 
phase convention for the helicity basis. In this paper, we adopt the quark helicity states 
and the 7-matrix representation of Bjorken and Drell ||. The anti-quark helicity states are 
chosen to be v(h, k) = Cu(h, k) T , where C 

differs from || convention). And thus, in the zero-mass limit, we have 
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Once again, the antiquark is moving in the —z direction. 

Clearly the coordinate system for treating quark polarization is just the one we have de- 
fined. On the other hand, the coordinate system associated with the anti-quark polarization 
(x, y, z) can be worked out from our choice for the antiquark helicity states and the definition 
for the axes in the last section. The z axis is clearly in the — z direction. Using Eq. (^), 
we find that the x axis is in the direction of the x axis and the y axis in the —y direction. 
[To avoid a kinematic zero, one has to restore the quark mass in the helicity states.] This 
coordinate system will be used for anti-particles throughout this paper. 

A straightforward calculation with the standard-model electroweak currents and the 
above helicity states yields, 
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where v e = 4 sin 2 6 W — 1 and a e 

electron to the Z. The couplings of the quarks to the Z are v u = 1 — |sin 2 ^^, Vd = v s 



+ | sin 2 9 W , a u 



1, and ad 



a s 
1 



Xi 



X2 



— 1 are the vector and axial vector couplings of the 
e quarks to the Z are v u = 1 — | sin 2 0w, v d = v s = 
— 1. N c — 3 is the color number of the quark, and 

(11) 
(12) 
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Eq. (|1~0|) indicates that the quark and antiquark from the same vertex have opposite helicity, 
a consequence of the massless limit. Therefore, as we will see in the next section, the polar- 
izations of A and A are correlated through the effects of spin- dependent quark fragmentation 
functions. 

The cross-section density matrix can also be expressed in terms of the Pauli matrices, 
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Any physical cross section can be obtained by taking the trace of the cross section density 
matrix with quark-antiquark polarization density matrices, which are defined as 



(14) 



For example, if we want to calculate the cross section when both quark and antiquark are 
polarized along the x axis, we have, 



da(e e + — > g(rr)c/(:r)) 



tr 



(icr 



(15) 



where .D 9 = F + cr 9 and .D 9 = F + cj are the density matrices for quark and antiquark 
polarization, respectively 



IV. FRAGMENTATION DENSITY MATRIX FORMALISM 

Quark fragmentation functions were introduced to describe hadron production from the 
underlying hard-parton processes. Apart from the well-known, spin-independent and chiral- 
even fragmentation function fi(z), there exist various chiral-odd and spin-dependent frag- 
mentation functions which are of particular interest because they describe novel spin effects 
in hadron production |2|||. At the leading twist, there are two additional fragmentation 
functions for spin- 1/2 baryon production, g%(z) and h\{z). All these fragmentation func- 
tions can be expressed as light-cone correlations in QCD, 

fi(z) = \zj ^e- a/2 (0|^(0)|A(P5)A) (A(PS)X$(\n)\0) , (16) 

9i{z) = \zj ^e-^ 2 (0|^(0)|A(PS,|)X) (A(PS {l )X\^Xn)\o) , (17) 

k(z) = \zj ^e-^(0| ? i75^(0)|A(P5 ± )x) (A(PS ± )Xm\n)\0) , (18) 
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where A(PS) represents the A hyperon with four- momentum P^ and polarization S^, 
normalized to P 2 = M A and S 2 = —M A respectively. We write P^ = p^ + n M M^/2, 

= Su + M A S^_, and Sjf = S ■ np^ + S ■ pn^ 1 with p^ and n M two null vectors (p 2 = n 2 = 
p ■ n = 1). The special component of p (n) is along (opposite to) the direction of the A 
momentum. These light-cone vectors will be defined separately for every observable hadron 
and will be labelled accordingly. The variable z, representing the momentum fraction of the 
quark carried by A, is defined in the usual way z = 2P ■ q/q 2 , where q is the momentum 
carried by the virtual boson. The summation over X is implicit and covers all possible 
states which can be populated by the quark fragmentation, and also the renormalization 
point (/i 2 ) dependence is suppressed. [QCD radiative corrections generate a different /i 2 de- 
pendence for each moment of these fragmentation functions, which is associated with the Q 2 
evolution of the experimental data. Although important for comparison with experiment, 
QCD evolution does not disturb the classification of spin dependent effects, so we suppress 
it throughout this analysis.] 

We turn to the density matrix formalism for quark fragmentation functions. Here we aim 
to construct quark and antiquark fragmentation density matrices that depend on helicity 
indices of quark (antiquark) and hyperon (antihyperon). Such a construction can proceed 
in three steps. First, we construct a density matrix in the Dirac representation. Then, we 
obtain a density matrix in the mixed representation, in which the helicity indices of hyperon 
(antihyperon) are manifest, but the quark (antiquark) is still labelled by the Dirac indices. 
In the final step, we transform the remaining Dirac indices into helicity indices. 

For the quark fragmentation q — > AX, we define 

M A {z,S,P) a p = J — e- lX / z (0\MV)\HPS)x)(A(PS)X\M\n)\0) . (19) 



Then, combining with Eqs. (|16|- |18|), we have 



M A (z,S,p/z) = ^t+^(n. S^+^j^J . (20) 

z z z 

to the leading twist. Another set of fragmentation functions can be defined to describe anti- 
quark fragmentation g -> AX. However, they can be related to the above fragmentation 
functions by using charge conjugation. In fact, if we define the anti-quark fragmentation 
density matrix as 

M- A {z,S,P) aP = |^e-^(0|^(0)|A(PS)x)(A( J PS)X|^(An)|0) , (21) 
we find by using charge conjugation, 

M A (z,S,p/z) = -C- 1 Ml(z,S,p/z)C , (22) 

where M A (z, S,p/z) means the transpose of matrix M A (z, S,p/z). Consequently for the 
ant i- quark fragmentation q —>■ AX, 

M- A (z,S,p/-z) = ^P~^(n.S^ 5 p + ^^J . (23) 



S 



Notice the sign change for the gi(z) term. 

The next step is to change M A , which depend functionally on the spin S* 1 of the A, to 
a 2 x 2 fragmentation density matrix carrying A (A) helicity indices. The density matrix is 
defined as, 

M(*,P)3f= J^e-^(0\M0)\HPH')x)(A(PH)X\MXn)\0) . (24) 

The transformation is possible because Eq. (|20|) with a general S contains all the information 
about the A spin dependence of the quark fragmentation. If we choose S to reproduce 
the helicity eigenstates, the diagonal elements of the spin density matrix are immediately 
obtained, 

M ++ = M A (z, = <, p/z) = + ^TSjf, 

Af__ = M A (z, = -e^p/z) = f -^i> - ^75^, (25) 

where we have chosen the direction of the A momentum as the z direction. [In the leading 
twist calculation, the direction of the quark and A can be taken to be collinear. The 
contributions from the A transverse momentum are among higher twists.] According to the 
superposition principle, we may extract the off-diagonal elements of the spin density matrix 
from Eq. (^). From 

|£" = M A e£> = ± (J) , |S"= M A e£> = -L Q , (26) 

we have 

M A (z, S» = M A e£,p/z) = ^(M ++ + M + _ + M_ + + M__) , (27) 

M A (z, S» = M A e^p/z) = \{M ++ + iM+- - iM-+ + M__) . (28) 

From these relations we can read off the off-diagonal components of M in the A helicity 
basis, 

M + - = ^75(^-^¥, (29) 

M- + = ^75(^ + ^¥. 

Here, the four-vectors and are defined by = (0,1,0,0), = (0,0,1,0) and 

= (P z , 0, 0, P°). In terms of the Pauli matrices {cx^} and 2x2 identity matrix E A , 
Eqs. (|25| ) and ( p9|) can be summarized as, 

M A (., p/z) = %^ E A + a A _ Mf) ^ + iy(J ^ . (30) 
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For the antiquark fragmentation, we define the density matrix analogous to Eq. (|24|), 



M(z,P)% s ' = |^e- lA / 2 "(0|^(0)|A(P^>(A(P^)X|^(An)|0) . (31) 
Then a totally parallel analysis shows that 

M- A (z, P/z) = ^ I A - ^4 of " ^75^ (44 + 4 A) > ( 32 ) 

A/ A, A/ 

with similar definitions 

M ++ = M A (z,^ = e^) , 

M__ = M A (^ = -e^) • (33) 

Here, the four-vectors e§, and e% are defined by = (0,1,0,0), = (0,0,1,0) and 
= (P z , 0, 0, P°) in the frame where the A momentum is in the z-direction. 
Finally we calculate the density matrices entirely in the helicity basis. This reformulation 
is straightforward — the Dirac matrix form of M\(z, pj z) is inserted between the helicity basis 
states u(h') and u(h). After taking care of the proper normalization, the resulting 2x2 
matrix, whose elements are labelled by the quark helicities, h'h, can be expanded in a quark 
Pauli matrix basis, as was just done for the A indices. The result has the remarkably simple 
direct product form, 

*!£±(z,p/z) = f x {z) F ® H A + 9l (z) a\ ® <x A + k(z) (o* ® a A + a« ® <x A ) . (34) 

It is clear from this form that §\ measures the probability that the longitudinal polarization 
of the quark transferred to that of A with fractional momentum z. The h\ measures the 
probability that the transversity of the quark is transferred. There is an obvious invariance 
of the hi term under rotations about the momentum direction, since transverse momenta 
have been integrated out of the fragmentation process. An analogous expression can be 
obtained for the anti-quark fragmentation into A, where the corresponding Pauli matrices 
will be defined in terms of the anti-particle's helicities. The direct product form for the 
antiquark case is 

^(z,p/z) = h(z) F ® I A + 9l (z) a! ® + h x (z) (of ® erf + a} ® af) . (35) 

The fact that Eqs. (0) and ( p5[ ) are identical is a direct consequence of charge conjugation 
symmetry, which is manifest in the formalism developed in §11. 

An expression like Eq. Q3~3] ) could have been written upon consideration of the rotation 
invariance and parity conservation of the parton fragmentation process. Aside from an 
explicit derivation, what is novel here is the identification of the coefficients with the leading- 
twist fragmentation functions, which have a formal connection with QCD. 
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V. A (A) DECAY DENSITY MATRIX 



It is well-known that the polarization of A and A can be measured through their weak 
decay. The relevant formalism and experimental information about their decay are included 
in the standard particle data book. Here, for our purpose, we need to obtain the spin density 
matrices for the decay. 

For the hyperon nonleptonic decay, the most general decay amplitude is 



M = G F m 2 u f (A - Bj 5 )ui 



(36) 



where A and B are constants and generally complex numbers, and Ui and Uf are Dirac 
spinors for the initial and final baryons, respectively. Analogous to Eq. (|36"D , the most 
general amplitude for the anti-hyperon nonleptonic decay is 



M' = -G F m 2 v t (A* + B* l5 )v } 



(37) 



with Vi and Vf Dirac spinors for the initial and final anti-baryons respectively. The spin 
density matrices for the A and A decay can be defined as, 



(38) 



D H 'h = "£{pir-\A{PH)){A{PH')\pir- 

Sp 

Ds>H = J2(pn + \A(PH))(A(PH , )\p7r- 



We have neglected the phase space integration which can be discussed separately. Substi- 
tuting Eq. (|3^) into Eq. (^) and again using the superposition principle, we find the four 
elements of D H > H , 



D ++ 
D + - 



(|a| 2 + |6| 2 )P P ■ P + (H 2 - \b\ 2 )m p M A + (a*b + ab*)P p ■ S R 
(H 2 + |6| 2 )P P • P + (H 2 - \b\ 2 )m p M A - (a*b + ab*)P p ■ S R 



2(a*b + ab*)M A P p ■ (e x - ie y ) 
2(a*b + ab*)M A P p ■ (e x + ie y ) 



(39) 



where Sr is related to the light-cone variables through Sr = p—n P p is the momentum 

of the proton, and 



a = G fm 2 + A 



b = -G f m3+B . 



(40) 



In terms of Pauli matrices, we obtain a compact form for the density matrix, 



D 



(\a\ 2 + \b\ 2 )P p -P + (\a\ 2 -\b\ 2 )m p M A 
+2(a*b + ab*)(P p ■ S R ) 
+2(a*b + ab*)M A \{P p ■ e x ) a£ + (P p ■ e y ) ai 



(41) 



A similar calculation for A gives 
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D 



\a\ 2 + \b\ 2 )P p -P + (\a\ 2 -\b\ 2 )m p M A 



-2(a*b + ab*)(P p - S R )a 



-2{a*b + ab*)M A (P p ■ e s ) a A + (P p ■ eg) a* 



(42) 



where P p is the momentum of the antiproton and Sr = p — fiM\j2. 

When performing the phase space integration, we keep the transverse momenta of the 
proton and the antiproton as differential variables. To take into account the particle decay 
width and the motion of the parent particle, it is convenient to define the boost-invariant 
quantities 



d 2 B 
d 2 B 



2M A r A 
1 



dP* 



d 3 P^ 



(2tt) 3 2E p (2tt) 3 2K 



dPl 



d z P, 



2M A r A J (2ir) 3 2E p (2n) 3 2E w 



D (2n) 4 S^(P-P p -P 7r ), 
D (2n)^\P-P p -P 7T ), 



(43) 



where T is the total width of the A decay and P n is the momentum of the pion. Particle 
masses have to be kept explicitly in the decay processes although they can be neglected in 
high-energy quark fragmentation. After integration, Eq. (H3T) becomes, 



d 2 B 



D 



d 2 P p ± 8(2vr) 2 M A r / 



p 

d 2 I) 



p-p+ 

v 

D 



P+P- 



d 2 P p L 8(27r) 2 M A r/ 



P-P p + - P+Pp 



(44) 



where P are defined in the usual light-cone coordinates, P ± = (P° ± P z )/\f2. 



VI. CROSS SECTION FOR e e+ -> qq -» AXAX -» pvr Xpvr+X AND ITS 

PHYSICAL SIGNIFICANCE 

Now we have the spin density matrices for all the subprocesses in e~e + — > qq —>■ 
AXAX -> pir~Xpir + X. Substituting Eqs. ([ID, ©, (H), and §D to Eq. (0), we get 
the final cross section 
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■a 



d£l K dzdz d 2 P±d 2 P± 



2(27r) 2 M A r/ 



x 



As 



p 

(\a\ 2 + \b\ 2 )(P p -P) ■ mm.M.s 



P + P- 



2(2vr) 2 M A r/ 



P-P+ _ p+p- 



q=u,d,: 



Q 2 (l + cos 2 0) + X2 (1 + cos 2 9)(t> 2 + a 2 q ){v 2 e + a 2 ) + Sv e a e v q a q cos 9 



-2Q 9 xi f e f 9 (l + cos 2 6) + 2a e a q cos 6 > ( fl{z)fl{z) + C gg gl{z)gf(z 



-i X2 



2v q a q (v 2 + a 2 )(l + cos 2 6) + 4t> e a e (t> 2 + a 2 ) cos 9 



-2Q qX i [a q v e (l + cos 2 9) + 2v q a e cos 9] } (C fg fl{z)gl{z) + C gf gl{z)ft 
Ql + X2(v 2 + a 2 e )(v 2 - a 2 ) - 2Q q xiv e v q ] cos(<p + tp) 
+2Q q Xi s Z _ ^ 2 v e a q sm(ip + <^)| sin 2 9 C hh h\(z)h\(z) 



(45) 



with ip and (p the azimuthal angles of the proton and the antiproton, respectively, in the 
coordinate system that is defined in the beginning of §111: the z axis is chosen to be the 
direction of the quark jet (or A) and the x axis in the plane of the beams and jets. The C's 
are given by 



a 



99 



a' 



: (P P ■ Sn)(Pp ■ Sr) 



Pn II" Pi 



M\ml 



a 



2 r p 



P\W 



mi 



r/r] 



C fg = a 



(Pp-S R ) 
M\m p 



\P 



a- 



m T , 



P\\ - 
V , 



(46) 



C, 



gf 



a- 



( Pp ■ Sr) _ \P p 

M A m p 



a- 



m r , 



Chh = of 



mt 



where P p ■ P = (Ml + m 2 — m 2 )/2 in C// and r\ {rf) is ±1 depending on whether the 
momenta of proton (antiproton) and A (A) are parallel or antiparallel. P p \\ and Pp» are the 
projections of the proton and antiproton momenta in the directions of A and A, respectively, 
in the respective rest frames of the parent particles, and P P L and Pp are projections of the 
proton and antiproton momenta onto the x-y plane. The a is the standard hyperon-decay 
parameter defined in the particle date table ]10| . 

There are three distinct classes of terms in the cross section. They correspond to three 
different type of angular dependences and are sensitive to different combinations of frag- 
mentation functions. The first class involves terms with products of fragmentation func- 
tions fi(z)fi(z) and gi(z)gi{z). They have no azimuthal dependence. The dependence on 
the polar angle arises entirely from the two-jet production cross section. To isolated the 
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gi(z)gi(z) term, one has to measure the correlation of the proton and antiproton momenta 
with respect to the momenta of A and A. The second class of terms contains the products 
of fragmentation functions fi(z) and gi(z); they arise from parity violation in electron or 
quark coupling with Z. The novel spin effects appear in the third class of terms, which 
involves the product of the transversity distributions h\{z)h\(z). Let us discuss the physics 
associated with this term in some detail. 

The origin of the sin 2 cos(<£> + (p) h\(z)h\{z) term is seen in the expression for the cross 
section into qq, Eq. (|i~3|), in the term proportional to {a% ® erf + (g) cr|). To appreciate the 
structure of that particular form, consider the annihilation through the pure photon channel. 
The intermediate state photon is produced with helicity of ±1, with equal probability for 
unpolarized incident leptons. If the qq pair were massive quarks produced near threshold, 
then their spins would be aligned with the photon spin direction, z.e.the beam direction. At 
production angle G near n/2, the term being considered must be maximum in magnitude. 
That requires the q and q spins to be preferentially parallel to one another and to align with 
the e~ or e + beam directions. Recall that the q and q momenta are in the ±£ direction and 
the x axis is in the scattering plane. Given the overall sign of the Q q term, it is seen that 
the cross section favors the o q x <g> erf expectation value be positive. Since the anti-particle's 
x-axis is oriented parallel to the particle's x-axis, this means that the x-components of q and 
q spin tend to be aligned. The corresponding y-components tend to be anti-aligned, since 
the anti-quark's f/-axis is antiparallel to the quark's y-axis in our convention. 

If the quark spin has an azimuthal orientation specified by <p q and the anti-quark, ipq, then 
the alignment of spins just specified leads to cos(0 g + <fiq) being positive, which favors the 
argument near 0. This result for the transverse components of the qq pair, can be visualized 
by having the q spin vector reflected through the scattering plane. Then the transverse spin 
vector of the quark tends to be parallel to the transverse spin vector of the reflected q, i.e. the 
quark transverse spin tends to line up with the mirror reflected anti-quark transverse spin. 

Next suppose that the quark (anti-quark) spin orientation is passed on to the A (A) 
fragmentation product. The asymmetry of the A decay into p7t~ provides a measure of the 
spin orientation. As Eq. (I"ID shows, the proton momentum tends to be aligned with the 
A spin. The corresponding Eq. fl42|) for the A decay into an antiproton yields the opposite 
distribution — the antiproton momentum tends to be antiparallel to the A spin orientation. 
Hence, while the x-components (y-components) of the quark-antiquark spins tend to align 
(anti-align), the x-components of the decay proton and antiproton momenta tend to anti- 
align (align). This is the interpretation of the cos((p + (p) term in the cross section of Eq. ([45]). 

Note that for the Z° intermediate state the conclusion is opposite from the above dis- 
cussion, because v q — a 2 is negative. So the transverse momentum of the proton tends to be 
aligned with that of the reflected antiproton. 

This construction leads to a simple phenomenological procedure for determining the value 
(for fixed z and z) of the product hx(z)hi(z). For the photon case, the above discussion is 
summarized by the statement that and Pp~ tend to be on the same side of the scattering 
plane. For the Z° case the tendency is for opposite sides of the scattering plane. So it is 
natural to define an asymmetry (for fixed z and z) via the number of proton-antiproton 
pairs on the same side of the scattering plane minus the number on opposite sides of the 
plane. The asymmetry selects the desired term, and has a simple B dependence, so that 
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results from all can be combined. The precise expression for this asymmetry follows from 
Eq. ©. 



VII. CONCLUSION 

In this paper, we have obtained a differential cross section for the process e~e + — > 
qq — > AXAX — > pir~Xpir + X. The cross section depends on, and therefore allows us to 
extract, the three twist-two quark fragmentation functions to the spin- 1/2 hyperon, fi(z), 
gi(z), and h\{z). These fragmentation functions contain important information about the 
soft QCD physics in quark fragmentation. Particularly interesting are the spin-dependent 
fragmentation functions which encode the behavior of the spin transfer. The cross section 
was obtained through a 2 x 2 spin density formalism in the helicity basis, in which the 
physics is made crystal clear. The formalism is general and can be used for other similar 
spin processes. 

The physical process that we discussed is accessible currently at LEP and SLAC. While 
up and down quark jets can produce A and A abundantly, the polarized hyperons are mostly 
produced from ss jets, as we expect from the constituent quark model. The production rate 
for ss jets is the same as dd and bb. The spin transfer from polarized s to A is expected 
to be large, especially in the large z region. The weak decay of A (A) into a proton and 
charged pion can be reconstructed easily and provides an excellent polarization analyzer. 
To obtain h\{z)h\{z) with enhanced statistics, an asymmetry can be defined by summing 
over all events with different Q and Pj_'s in bins of z and z. Given these remarks, we are 
looking forward to a first measurement of the spin-dependent fragmentation functions! 
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FIGURES 

FIG. 1. Feynman diagram for the process e~e + — ► qq — > AX AX — > pir~ Xpir + X . The arrows 
denote directions of momentum flow (the fermion number flow for antiparticles is against the 
momentum flow). The helicity indices are explained in the text. 
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